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ABSTRACT 

We exhibit soliton solutions of QCD in two dimensions that have the quantum 
numbers of quarks. They exist only for quarks heavier than the dimensional 
gauge coupling, and have infinite energy, corresponding to the presence of a string 
carrying the non-singlet color flux off to spatial infinity. The quark solitons also 
disappear at finite temperature, as the temperature-dependent effective quark 
mass is reduced in the approach to the quark/hadron phase transition. 
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1. Introduction 



One of the key problems in non-perturbative QCD is the understanding of 
"constituent" quarks from first principles. Historically, the first indications for 
the physical reality of quarks came from the successes of the constituent quark 
model in light hadron spectroscopy and matrix elements, assuming m U} d — 300 
MeV and m s ~ 500 MeV. 111 Subsequently, current algebra was abstracted from 
the lagrangian for the light quarks (u, <i, s) and applied successfully to calculate 
the properties and interactions of pions and kaons. However, the success of 
quark current algebra could be understood only in the context of approximate 
chiral symmetry, according to which the u, d and s masses in the lagrangian were 

[3] 

much smaller that the original "constituent" masses. The so-called "current 
quark" masses are a few MeV for u and d and (9(100 to 200) MeV for the s 
quark.' 41 Thus the question arises: what is the relationship between "current" 
and "constituent" quarks in QCD? In the cases of the "heavy" quarks (c, 6, t), 
the difference between "constituent" and "current" quarks is presumably not so 
great, although even the successful quarkonium potential model has not been 

[5 — 9] 

derived completely from QCD. Even more of a challenge is to derive the 
constituent versions of the light u, d and s quarks from QCD. 

Refs. [6], [8] proposed that the constituent quarks are soliton solutions of a 
hypothetical effective lagrangian of QCD, with effective bosonic fields carrying 
both color and flavor. It is not known whether such an effective lagrangian exists 
in a well defined sense, nor is it known whether there are stable, non-trivial 
classical solutions. Sidestepping these questions, ref. [8] makes the interesting 
observation that the conjectured group-theoretical structure favors solitons with 
quantum numbers of quarks. In the absence of a well-justified effective dynamics 
of QCD, however, further progress requires a theoretical laboratory where the 
relevant phenomena are explicitly calculable. 

A useful laboratory for studying this problem is QCD in 2 dimensions.' 10 201 
This theory can be written in bosonized form' 141 for arbitrary numbers of colors 
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N c and flavors Nf. It reflects accurately the phenomena of quark confinement 
and condensation in the vacuum that we expect to occur in QCD in 4 dimen- 
sions. Even though Goldstone bosons cannot exist in 2 dimensions, QCD 2 yields 
a relation between the pion mass, light quark mass and quark condensate sim- 
ilar to that expected on the basis of approximate chiral symmetry in QCD 4 . [111 
Moreover, QCD 2 has finite-energy soliton solutions for arbitrary values of N c and 
Nf that can be interpreted as baryons) 16 191 in close analogy with the Skyrmion 
interpretation of baryons as solitons in QCD 4 . [21 231 

In this paper we start an investigation of constituent quarks in QCD 2 . Specif- 
ically, we examine QCD 2 with one heavy quark flavor, and look for soliton solu- 
tions corresponding to a single quark, i.e. baryon number 1 (the baryon number 
is normalized to be N c for the nucleon). They exist, for a sufficiently heavy quark 
<5, but have infinite energy, corresponding to a string carrying the non-singlet 
color flux off to spatial infinity. These quark soliton solutions disappear when the 
meson mass parameter M is reduced to become comparable to the gauge cou- 
pling strength e c (which, we recall, has the dimension of mass in QCD 2 ).They 
disappear when it is energetically favoured to create a new QQ pair, rather than 
to create a string with tension (a/) -1 ~ e 2 c (for the analogous result in QED 2 , see 
ref. [24]). The quark solitons also disappear when the temperature T is increased 
to O(M), because of a reduction in the T-dependent effective meson mass. The 
quark condensate and the baryonic solitons disappear at a somewhat higher tem- 
perature, signalling deconfinement at the quark/hadron phase transition. 
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2. Quarks as Solitons 



We start with the bosonized effective action for QCD 2 



[16,17] 



A(N F =1) = S[h] + ±J d 2 x{d^f 



+ M / d 2 x Tr 



e 

2tt 



/ d 2 xTrH 2 



(1) 



where we took one flavor Nf=l for simplicity, N c is the number of colors, h is 
the SU(iV c ) color matrix of field operators, and 

S[h] = ^ J d 2 xTr (d^hd^h)) + -^-J d 3 ye ljk Tr (h)d l h) (h^djh) (h)d k h) 

(2) 

where the last term is the Wess-Zumino term. In eq. (1) e c is the dimensional 
gauge coupling and <fi is the baryon number field. The mass scale M is related 

to the current quark mass tuq and to the gauge coupling, in the limit where 

^ u [16] 
e c > m Q , by 



M = [CN c m Q (^=y] — . p = i _ 
V 71 " 



(3a) 



where C 



\e 1 ~ 0.89, and 7 is Euler's constant. We will be interested in the 



opposite limit tuq ^> e c , hence we choose M—tuq. The field H in (1) is related 
to h through 

d-H = ihd-h + (3b) 

We define hermitian fields ip by 



h = e 1 ^* 



(4) 



and take Tr</? = 0, absorbing a possible trace term Tr</? = , n in 



4tt 
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To search for the lowest-energy classical solutions,' 161 we take (p diagonal with 
entries (pk, (k = 1, 2, ... , N c ). 



k=i 

+ 2M 2 J2 j cos (yfotpk + J^-A - 2e 2 c J2 f 
k=i J ^ V c / k=1 J 



Defining the shifted variables 



Xk = <Pk + 



{k = l,2,...,N c ) 



we find that for static Xk the Hamiltonian takes the form 



H (Xk) = \Y.j (diXk) 2 + V 

k- 1 

V = 2M 2 



k=i 



1 — COS V 4:TCXk 



fc=i 



2 2e| 

Nr. 



J2xi 
i=i 



and the baryon number is' 16 ' 171 



1 Nc 

B = —j= V[xfc(oo) - Xk(-oo)] 



The field equations for the static case are 



Xk - 4M 2 ^ sin V4nxk - 4e 



Xk 



--Y 



Xi 



1=1 



(5) 



(6) 



(7) 



(8) 



= (k=l,2,...,N c ) (9) 



Let us look for solutions where Xk(— °o) = for all k. Then, at x = oo, we must 
have 



M 2 v / 7r sin v / 47rx/ £ (cx)) + e 



1 ^ c 

Xk(oo) -jrJ2xi(oo) 
c i=i 







(10) 



If we seek a solution with Xfc(oo) = %(oo) for all k, we find that Xfc(°°) = 
^A/7rn, and £? = \nN c . But to have the eigenvalues of the squared mass-matrix 
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[d 2 V/dxkdxi] all positive, we must have even n, and thus integer baryon number 
B=N C , corresponding to a nucleon state. 

In our search for quark solitons, let us first consider the simplest case N c = 2. 
Then, the equations (10) become 

sin v^X^oo) = -e v / 7r(xi(oo) - % 2 (oo)) (11a) 

sin v / 4ti : X2(oo) = -ev^ (% 2 (oo) - xi(oo)) (lib) 

where e = e^/27rM 2 . Looking at the equations (11a), (lib), we see that we can 
have 

X 2 (oo) = -Xi(oo) + ti-s/tt (12a) 



or 



X 2 (oo) =xi(oo) + (n + ^)v^ (12b) 



i 

In case (a) we have B=n, and in case (b) we have B = -j=xi(°°) + ( n + 2 )- 
However, looking again at the eigenvalues of the matrix of second derivatives 
of the potential, and requiring positivity, we exclude case (b) (which otherwise 
would have non-integer baryon number). 

For case (a), where B = n, we may have quarks for n = 1, and other non- 
baryonic solitons for odd values of n. As seen in the figure, there is no such 
solution when e = e 2 c /2^M 2 > 1 (for the analogous result in QED 2 , see ref. [24]). 
This can also be seen directly from the fact that |sin(x)/x| < 1 for any x. For 
e « 1, however, we find a series of solutions with positive second derivative 
matrix* Defining 



Xi(oo) - \ n V^ 



(13) 



* See ref. [26] for comparison. 
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these are 

(2/) for n even 



6 = 



e 2 

7T 



2M 2 



(21 + 1) for n odd 



(14) 



2 

in the limit where Ij^i <^ ^" ^ ne number of solutions increases as e = e 2 /2nM 2 
decreases, as seen in the figure, and one can find solutions numerically for the 
general case e < 1. The solutions (14) above have infinite energy, with classical 
string tension 

„ f (2/) 2 for n even 
(a')- 1 ^ 2 ^ (15) 
\ (2/ + l) 2 for n odd 

They correspond to excitations of "colored" states. The case n=l, (2/+l)=l is 
the single constituent quark soliton. 

Note that the string tension (15) grows like I 2 for even n, and like (/ + ^) 2 for 
odd n. A possible interpretation is that the classical string tension is growing like 
J 2 , where / is the color isospin of the state: quantum mechanically it obviously 
has to become /(/ + 1). 



3. Generalization to N c > 2 

It is a simple matter to generalize the above asymptotic solution (12a) to an 
arbitrary number of colors. The equations (10) clearly require 

\fAnXk{°c) = n k^ + O { e c/ V^M 2 ) for some integers {n^}. The positivity of the 
eigenvalues of [d 2 V/dxkdxi] in fact requires that the {nk} be even: nk = 2p k , 
where the {pk} are integers. 

e 2 

Hence we look for asymptotic solutions, for ^—^ <C 1, of the form 

2ttM 2 

V^x*(oo) = 2p k n + Ck ^ + (J^J (k=l,...,N c ) (16) 
The equations (10) in fact tell us that for integer {p k } the corrections {c^} are 
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given by 



* (§»)) 



c k = 2\-p k + — | y jPj \ | (fc = l,...,JV c ) (17) 



Equation (8) then tells us that the baryon number 



i v ' v 



-^Ex*(oo) = X)p* (18) 

V k=l k=l 



as ^Cfc = 0. Hence B is an integer. The variable £ (13) for the case iV c =2 may 



k=i 

be generalized to the variables 



which take the values 



2/ 1 ^ \ / 2 \ 2 

e fc = (2p*-s)7r-^u-— X; Pl J + o^j (fc = i,...,jv c ) (20) 

generalizing the expression (14) for N c =2. Correspondingly, the string tension 
(15) becomes 

(a')- 1 =2ne 2 c J2(B/ N c-Pk) 2 (21) 



fc=i 



We see that when £ is some integer multiple of N c , corresponding to a multiple- 
baryon state) 181 the string tension vanishes, (a/) -1 = 0, if all the {pk} are taken 
equal to the multiple-baryon number. 
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4. Finite temperature 



Further insight into the interpretation of these QCD2 solitons can be gained 
by considering the bosonized form of the action (1),(5) at finite temperature. The 
temperature-dependent one-loop corrections to a generic 2-dimensional action 
can be written in the form 120 ' 



00 

SV (T) = L I dk Tr In - e ~^^ + M^ (22) 




where T = 1//3 is the temperature and M is the mass matrix. The leading 
field-dependent term in (22) has the form 

SVi (T) = c T Tr (23) 

where c is a coefficient independent of the specific theory. In the case N c = 2, the 
finite-temperature correction (23) to the effective potential (7) takes the form 

SVi(T) = ct(^/2m(^[cosV^xi + \J cos v / 4^% 2 ^ 

el ( 1 1 " < 24 > 

+ = + 



&\[2t; V Ma/cos v^Xi Ma/cos v^X 



[25] 

We now discuss some likely implications of this correction whilst recognizing 
that a complete discussion of the phase structure of QCD 2 would require evalu- 
ating the effective potential to all orders in the loop expansion. 

As discussed in Section 2, "baryonic" soliton solutions to the zero-temperature 
N c = 2 field equations (even B) exist for arbitrary values of e = e 2 c /2nM 2 , whilst 
"quark" (odd B) solitons exist only for e < 1. As follows from equations (13), (14), 
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they correspond to 



Xi(oo) ~ m^/ir + e^^ -mj (m e J) 



(25) 



in the limit e < 1, and ^2(00) = — Xi(oo) + n^pK. In the case of even B, there 
are baryonic solitons with m = B/2 for which the e-dependent correction in (25) 
vanishes, whereas it is always present for quark solitons with odd B, and indeed 
becomes large when e — > 1, corresponding to the previously noted disappearance 
of the soliton. When e C 1, the leading effect of the finite-temperature correction 
(24) is to replace the zero-temperature solution (25) by 



Xi (00) ~ mVir + e(T)V^Tv ^ 



f-m) (20) 



where 



£(T) i 1 - 571 ) (27) 

The effect of the finite-temperature modification (27) is in the direction of desta- 
bilizing the quark solitons (and also baryonic solitons with m 7^ B/2), which 
presumably disappear at critical temperatures T cB m = O(M), with solitons 
whose values of \B/2 — m\ are largest disappearing first. In particular, the last to 
disappear would be the quark solitons with \B/2—m\ = |, and finally the baryons 
with B = m. The latter actually disappear only when the coefficient of the peri- 
odic contribution to the zero-temperature potential (7) is finally overwhelmed by 
the finite-temperature corrections (24), which occurs at T ch B j 2 = T q / h = O(M). 

We now discuss briefly some possible interpretations of these results. We in- 
terpret Tq/h ~ 0(M) as the quark/hadron phase transition temperature, where 
the fermion condensates dissolve,' 27 ' 281 , the baryonic solitons disappear, and cor- 
respondingly quarks are deconfined, as discussed in refs. [29], [30]. We are not 
surprised by the fact that T q / h seems to be quark mass- (and hence flavor-) de- 
pendent in QCD 2 , whereas it is expected to be O(Aqqd) and flavor-independent 
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in QCD4. In QCD4, Tg/h can be determined by entropy considerations (the flux 
tube fluctuates more and more until flux lines eventually fill all space) , whereas in 
QCD2 the flux line cannot fluctuate out of the single space dimension, and T q / h 
is determined by the energetics of quark pair-creation, which becomes favoured 
when T = 0{M). The presence of fermion condensation is necessary for solitons 
to form, but not sufficient . Even when T = 0, for any fixed value of e solitons 
with large values of \^ — m\ do not exist and the number of such solitons dimin- 
ishes as T increases, with the "quark" solitons disappearing before the B = 2m 
"baryons" . 

We plan to return to these and other issues at finite temperature and chemical 

[31] 

potential in a future publication 

5. Conclusions 

We have shown in this paper that QCD2 has quark soliton solutions if the 
quark mass is sufficiently large. We have also shown how these quark solitons 
disappear when the quark mass tuq is reduced until the meson mass M (3a) 
becomes comparable to the dimensional gauge coupling strength e c , or when the 
temperature is increased to 0{M) in the approach to the quark/hadron phase 
transition. The next step in this approach to the derivation of constituent quarks 
from QCD is to look for meson solutions to the QCD 2 field equations which 
contain one heavy and one light quark. The light meson cloud in the presence of 

a heavy quark would correspond to the concept of a light constituent quark in 

[32] 

QCD. Once these solutions are obtained, we hope to be able to abstract the 
relevant features of the field equations in QCD4 and then solve them to construct 
constituent quark solitons also in 4 dimensions. 
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FIGURE CAPTION 

Comparison of the left- and right-hand sides of the soliton eq. (11a), applying the 
relation (12a) corresponding to a quark soliton with n = 1, for e = e 2 ./2nM 2 = 0.1 
(solid line) and e = 0.5 (dotted line). There are no solutions for e > 1 (dot-dashed 
line, drawn for e = 2). 
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